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EXECUTIVE  SUMMARY 


OBJECTIVE 

The  objective  of  this  research  was  to  explore  the  effects  of  altering  the  third  junction  in  an  array  of 
bi-superconducting  quantum  interference  devices  (bi-SQUIDs)  on  the  overall  performance  characteris¬ 
tics.  The  third  junction  could  be  altered  by  either  fabricating  it  using  a  different  method,  such  as  ion  beam 
milling  or  ion  damage,  or  with  junctions  specifically  designed  to  exploit  certain  properties,  as  with  multi- 
ferroic  junctions.  First,  an  exploration  was  performed  to  review  the  effects  of  capacitance  and  noise  terms 
for  a  single  junction,  and  a  single  SQUID.  This  was  done  to  validate  the  simulations  and  to  create  a  base¬ 
line  with  which  to  compare  the  new  work.  Once  the  effects  were  characterized,  a  single  bi-SQUID  and 
array  of  bi-SQUIDs  were  examined  that  had  noise  temperature  and  capacitance  terms  on  the  bisecting 
junctions  only.  This  setup  represents  a  device  with  traditional  main  junctions  and  different  third  junction, 
which  could  be  used  to  precisely  design  the  array  for  tuneabilty  or  filter- like  effects  by  controlling  the  lin¬ 
earity  and  voltage  swing  of  the  anti-peak.  The  linearity  and  voltage  swing  are  important  for  increased  gain 
and  improved  accuracy.  Finally,  a  frequency  dependence  was  added  to  the  third  junction  to  determine  the 
effects  of  having  a  multiferroic  third  junction. 

RESULTS 

We  determined  that  for  a  single  junction  and  single  SQUID,  capacitance  is  detrimental,  as  it  causes 
hysteresis.  Even  a  small  amount  of  hysteresis  in  the  I-V  curve  destroyed  the  voltage  swing  in  the  V (xg) 
curve,  which  is  the  response  we  are  trying  to  manipulate  for  signal  detection.  A  decrease  in  the  junction 
resistance  counteracts  the  hysteresis  arising  from  the  capacitance.  A  large  resistance  negatively  impacts 
junction  performance  and  has  an  even  stronger  impact  on  hysteresis  since  f3c  oc  while  /3c  oc  C.  A  small 
amount  of  temperature  noise  also  decreases  any  hysteresis  present;  however,  it  also  decreases  the  junction 
critical  current.  When  hysteresis  is  present  in  the  I-V  curve,  the  bias  current  can  be  decreased  to  increase 
the  voltage  swing;  however,  this  results  in  hysteresis  in  the  V (xg)  curve. 

When  we  just  include  the  capacitance  and  temperature  terms  on  the  bisecting  Josephson  junction, 
and  leave  the  outer  loop  main  junctions  as  ideal  junctions,  their  impact  on  the  V (xg)  curve  are  decreased. 
The  temperature  term  counteracted  the  effects  of  a  large  third-junction  critical  current  on  the  V  (xg)  curve 
when  there  were  small  capacitance  and  resistance  values.  Gain  is  increased  as  we  couple  the  SQUID  or  bi- 
SQUID  devices  together  into  arrays,  so  when  we  fabricate  the  prototype  bi-SQUID  sensor,  we  anticipate 
having  tens  of  thousands  of  bi-SQUIDs  in  the  array.  As  with  the  single  bi-SQUID,  the  impact  of  the  ca¬ 
pacitance  and  temperature  terms  on  the  V  (xg)  curve  was  decreased.  From  this  we  can  surmise  that  if  care 
is  taken  with  the  main  two  Josephson  junctions,  then  the  bisecting  junction  can  be  fabricated  with  an  infe¬ 
rior  method  and  still  result  in  a  good  device.  An  inferior  method  could  be  used  to  overcome  design  layout 
issues  or  for  the  properties  gained  by  using  them  (i.e.,  multiferroics). 

For  the  frequency  dependence,  we  determined  that  we  can  control  the  signal  detection  if  there  is  a  fre¬ 
quency  dependence  on  the  critical  current  of  the  third  junction.  For  a  very  basic  example,  if  we  could  find 
a  method  that  results  in  a  large  ics  value  10)  for  large  frequencies,  and  have  a  value  close  to  ics  =  1.0 
in  a  range  where  the  signal  strength  is  small,  and  ics  =  0  elsewhere,  we  could  have  a  device  that  acts  like 
a  low-pass  filter  and  amplifies  the  signals  further  in  the  determined  range.  The  value  of  the  third-junction 
critical  current  had  a  much  larger  effect  on  the  frequency  dependence  than  the  changes  to  the  time  con¬ 
stant. 
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RECOMMENDATIONS 


Based  on  the  results  of  this  study,  we  determined  that  a  modification  to  the  third  junction  is  a  viable 
approach  to  introduce  a  desired  frequency  dependence  and/or  control  of  the  linearity  and  voltage  swing  of 
the  anti-peak.  This  study,  while  extensive,  was  by  no  means  exhaustive,  and  a  significant  amount  of  work 
remains  to  fully  examine  the  different  effects. 
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1.  INTRODUCTION 


Superconductive  materials  have  the  unique  property  that  they  lose  all  resistance  to  electric  current 
when  cooled  below  a  critical  temperature.  A  Josephson  junction  is  a  thin  layer  of  insulating  material  sep¬ 
arating  two  superconductors  that  is  thin  enough  for  electrons  to  tunnel  through.  Two  Josephson  junctions 
placed  in  parallel  into  a  ring  of  superconducting  material  form  what  is  known  as  a  DC  Superconducting 
Quantum  Interference  Device  (SQUID).  SQUIDs  can  sense  minute  magnetic  fields  approaching  10“^^ 
Tesla.  These  SQUIDs  can  be  arranged  in  arrays  with  different  coupling  schemes  and  parameter  values  to 
create  a  combined  voltage  flux  response  (V (xe)).  Control  over  the  linearity  of  V (xg)  is  important  for  in¬ 
creased  gain  and  improved  accuracy. 

Increasing  the  number  of  SQUIDs,  N,  in  an  array  increases  the  detected  signal  amplitude  as  a  function 
of  N,  while  the  noise  is  increased  as  a  function  of  [1].  A  single  SQUID  has  two  Josephson  junctions 
in  a  ring  of  superconducting  material.  A  third  junction  bisecting  the  superconducting  loop  is  added  to  in¬ 
crease  the  linearity  of  the  voltage  response  output  [2].  This  device  is  called  a  bi-SQUID.  We  will  explore 
using  the  effects  of  using  a  different  material  and/or  method  on  the  bisecting  Josephson  junction  for  high 
temperature  superconductor  (HTS)  YBa2Cu307_5  (YBCO)  bi-SQUIDs. 

This  technical  report  is  organized  as  follows.  First,  we  review  the  effects  of  capacitance  and  noise 
terms  for  a  single  junction  [3-6].  We  then  expand  the  analysis  to  a  single  SQUID  [7].  Next,  we  compare 
the  single  SQUID  with  a  single  bi-SQUID  that  has  noise  temperature  and  capacitance  terms  on  the  bi¬ 
secting  junction  only.  Then  we  explore  the  effects  of  noise  temperature  and  capacitance  on  an  array  of 
50  bi-SQUIDs.  Finally,  we  look  at  effects  of  frequency  dependence  of  the  third  junction.  The  effects  ex¬ 
plored  are  those  that  could  arise  from  junctions  either  fabricated  using  a  different  method,  such  as  ion 
beam  milling  or  ion  damage  [8-11],  or  with  junctions  specifically  designed  to  exploit  certain  properties 
as  with  multiferroic  junctions  [12-14].  Precisely  designing  the  junctions  could  result  in  benefits  such  as 
tuneabilty  or  filter-like  effects. 
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2.  SINGLE  JUNCTION  AND  HYSTERESIS 

First,  we  examine  the  single  Josephson  junction  dynamics  and  the  effect  of  capacitance,  resistance, 
and  temperature. 

2.1  MODELING  EQUATIONS 

In  Figure  1  we  show  the  schematic  diagram  of  a  single  Josephson  junction.  We  model  a  single  resis- 
tively  and  capacitively  shunted  Josephson  junction,  where  the  X  represents  a  Josephson  junction,  Ii  is  the 
current  through  the  junction,  and  /b  is  the  DC  bias  current. 

/. 


(p 

Figure  1.  Single  Josephson  junction  schematic. 


The  equation  to  model  the  voltage-current  response,  which  includes  a  temperature  noise  term,  is 


nc 

2eIo 


"I  ft  .  h  Ic  .  , 


where  Ti  is  Plank’s  constant,  C  is  the  junction  capacitance,  e  is  the  charge  of  an  electron,  Iq  is  the  junction 
critical  current,  R  is  the  junction  resistance,  bb  is  Boltzmann’s  constant,  T  is  the  temperature  in  Kelvin, 
Ic  is  the  junction  critical  current  accounting  for  fabrication  imperfections,  (f  is  the  phase  difference  across 
the  junction,  and  5{t  —  t')  is  a  Gaussian  noise  delta  function.  Other  terms  we  are  disregarding  for  the  mo¬ 
ment  are  the  second  harmonic  (q  sin  2(p)  and  the  AC  drive  current  (Id  sin  ujdt)  [6].  The  AC  drive  current 
results  in  a  phenomenon  known  as  Shappiro  steps  [15].  For  our  application,  we  only  apply  a  DC  drive  cur¬ 
rent.  Setting  T  =  Josephson  time  constant)  gives  the  dimensionless  ordinary  differen- 

Ti 

tial  equation  (ODE): 


dip 


.  .  ,  ^ 

icSmcp  H - —— — Zd[T  - 

IqTi 

2eC  R^I 

where  /3c  = - ^ - is  called  the  McCumber  parameter  and  is  found  by 


hc  2  _  /2eRIoy  _  2eCR‘^Io 

2^^°  ~  2^  V  ^  / 


h 

From  the  drawing  and  Kirchhoff’s  current  law,  we  have  that  ii  —  is,  so  the  second  order  differential  equa¬ 
tion  we  need  to  solve  to  get  the  current  voltage  (I-V)  curve  is 

_  d?ip  dip  ,  .  .  f. 


where 


n  .  dp>  dp> 

1/  =  -V  = 


Splitting  the  second  order  differential  equation  into  two  dependent  first-order  differential  equations  yields 


=  V 


dip 
dr 

n  dv  .  .  .  f. 

PC-r~  —  ^B“'^cSin(^H - - — 2o(t  —  t)  —  v. 


dr 


hn 


(1) 
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Equation  (1)  is  integrated  in  Matlab  using  the  Euler-Maruyama  algorithm  with  a  time  span  of  r  =  (0, Ar 
over  a  range  of  values  from  negative  to  positive  and  then  sweeping  back  positive  to  negative.  This  time- 
dependent  voltage  response  output  is  then  filtered  through  a  moving  average  filter  and  plotted  against  is- 


2.2  PARAMETERS 


The  parameters  related  to  the  modeling  of  Equation  (1)  are  summarized  below: 
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2ph  =  6.6260696  x  10“^^ Js 

^  =  2.067834  X  10“^*^Wb 
2e 


Rio 

2ei?/o 


n 

2eCR^lQ 

n 

T 

—  X  10 
^0 


t  =  f  =  27rVQ^Qt 


$0 
27iCR^h 
$0 


2^{RC){hR) 

$0 


-7 


2.3  SINGLE  JUNCTION  SIMULATION  RESULTS 

A  small  amount  of  temperature  noise  (also  known  as  thermal  noise)  and  capacitance  with  a  normalized 
amount  of  resistance,  in  a  single  junction  produces  the  behaviors  as  shown  in  Figure  2.  Here,  Kp  —  0.01, 
i?  =  1.0,  and  C  =  0.1.  The  plot  on  the  left  is  the  voltage  over  time  of  the  junction  as  the  bias  current  ip 
is  varied  from  —4.0  to  4.0  and  back  to  —4.0.  The  plot  on  the  right  is  the  voltage  plotted  against  ip  using  a 
moving  average  filter.  There  are  no  observable  effects  of  noise,  and  no  hysteresis.  This  is  an  example  of  an 
ideal  junction’s  dynamics. 


Figure  2.  V{t)  curve  (left)  and  I-V  curve  (right)  with  Kp  =  0.01,  R  =  1.0,  and  C  =  0.1. 

Increasing  the  capacitance  C  increases  the  amount  of  hysteresis  present  in  the  I-V  curve.  For  Kp  — 
0.01,  R  —  1.0,  and  C  —  2.0  (see  Figure  3).  The  amplitude  is  decreased  in  the  time  series  plot  of  the  volt¬ 
age  on  the  left.  The  hysteresis  in  the  I-V  curve  on  the  right  allows  you  to  see  the  time  progression  color 
change.  At  time  zero,  the  line  is  red  then  goes  to  orange,  yellow,  green,  and  ends  at  blue.  The  hysteresis 
occurs  when  the  current  increases  from  negative  to  positive  and  delays  when  the  response  enters  the  su¬ 
perconductive  state.  A  current  with  magnitude  smaller  than  the  critical  current  is  needed.  Once  there,  the 
junction  stays  in  at  zero  voltage  until  the  positive  critical  current  is  reached.  On  the  path  from  positive  to 
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negative,  the  same  thing  happens  but  this  time  the  delay  is  on  the  positive  side  and  the  junction  stays  at 
zero  voltage  until  the  negative  critical  current  is  reached. 


Figure  3.  V{t)  curve  (left)  and  I-V  curve  (right)  with  Kp  =  0.01,  R  =  1.0,  and  C  =  2.0. 


When  the  capacitance  is  increased  further,  the  oscillation  amplitude  is  decreased  further  to  almost  non¬ 
existent  and  the  hysteresis  is  increased,  and  shown  in  Figure  4  where  Kp  =  0.01,  R  =  1.0,  and  C  =  20. 
The  hysteresis  is  an  undesirable  quality  in  the  I-V  for  the  applications  in  which  we  are  interested.  As  such, 
it  can  be  tentatively  concluded  that  a  high  capacitance  is  bad. 


Figure  4.  V{t)  curve  (left)  and  I-V  curve  (right)  with  Kp  =  0.01,  R  =  1.0,  and  C  =  20. 


Changing  values  of  R  changes  r  and  fic  since  both  parameters  depend  on  R,  where  r  (x  R  and  fic  ^ 
R^.  Decreasing  R  below  one  when  /Ic  is  small  makes  the  solution  unstable  for  sufficiently  small  values  of 
R.  This  is  because  the  inverse  of  the  matrix  used  in  the  modeling  becomes  practically  singular  when  the 
value  of  fic  is  very  small.  The  phenomenon  is  due  to  modeling  instabilities  rather  than  what  occurs  in  a 
fabricated  junction;  there  are  alternate  model  equations  for  when  jSc  —  0. 

For  larger  /Ic,  reducing  R  reduces  the  hysteresis;  see  Figure  5  where  Kp  =  0.01,  R  =  0.5,  C  =  2.0, 
and  flc  —  ^  (0.5)^  =  0.5.  For  these  values,  the  response  is  approaching  the  no  hysteric  ideal  solution. 

There  is  a  slight  dampening  of  the  amplitude  of  oscillations  for  large  negative  and  positive  values  of  the 
current  and  voltage.  So,  the  hysteric  effects  of  a  large  capacitance  can  be  mitigated,  if  it  is  accompanied  by 
a  small  R. 

When  non-negligible  capacitance  is  present,  increasing  the  value  of  R  increases  the  junction  hystere¬ 
sis.  The  hysteresis  is  increased  faster  than  when  the  capacitance  is  increased  due  to  R  being  squared  in 
/3c7.  Figure  6  shows  the  result  of  very  large  hysteresis,  where  Ffp  =  0.01,  i?  =  5.0,  C  =  2.0,  and 

=  2  X  (5.0)^  =  50.  When  the  current  is  negative  increasing  to  positive,  the  superconductive  state  is 
not  achieved  until  the  current  becomes  zero,  then  the  voltage  stays  zero  until  the  critical  current  is  reached. 
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From  positive  to  negative  current  again  the  response  doesn’t  go  into  the  superconductive  state  until  the  cur¬ 
rent  reaches  zero,  then  the  voltage  stays  zero  until  the  negative  critical  current  is  reached. 


Figure  5.  V{t)  curve  (left)  and  I-V  curve  (right)  with  Kp  =  0.01,  R  =  0.5,  and  C  =  2.0. 


Figure  6.  V{t)  curve  (left)  and  I-V  curve  (right)  with  Kp  =  0.01,  R  =  5.0,  and  C  =  2.0. 

For  the  Josephson  junction  with  =  1.0,  i?  =  1.0,  and  C  =  2.0  modeled  in  Figure  7,  the  param¬ 
eter  R  has  been  decreased  to  the  normalized  value  and  the  noise  temperature  has  been  increased.  There 
is  some  hysteresis  present  since  C  —  2.0,  but  a  reduced  amount  than  in  the  previous  case.  The  transition 
point  (critical  current)  and  hysteresis  are  both  decreased.  There  are  much  greater  amplitude  oscillations  for 
Kp  —  1.0,  especially  while  in  the  superconductive  state,  when  compared  with  all  the  other  examples  so 
far. 


Figure  7.  V{t)  curve  (left)  and  I-V  curve  (right)  with  Kp  =  1.0,  R  =  1.0,  and  C  =  2.0. 


If  we  plot  the  previous  example  where  Kp  —  1  with  the  example  from  above  with  the  same  parame¬ 
ters  except  with  Kp  =  0.01,  and  zoom  in  to  the  positive  current  hysteresis  loop  we  get  the  plots  in  Figure 
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8.  Black  is  the  warmer  temperature.  It  can  be  seen  clearly  that  while  the  hysteresis  decreases,  the  critical 
current  also  decreases. 


Finally,  for  much  larger  temperatures,  the  superconductive  state  disappears  completely.  Figure  9  has 
Kp  =  10.0,  R  =  1.0,  and  C  =  2.0.  Even  when  the  current  is  zero,  there  is  sufficient  temperature  noise  to 
keep  the  junction  in  the  normal  state.  Now  that  we  have  explored  the  dynamics  of  the  basic  building  block 
of  SQUID  arrays  we  will  move  on  to  a  single  SQUID,  which  is  a  ring  of  superconducting  materials  with 
two  Josephson  junctions  placed  symmetrically  across  from  each  other. 


Figure  9.  V{t)  curve  (left)  and  I-V  curve  (right)  with  Kp  =  1.0,  R  =  10,  and  C  =  2.0. 
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3.  SINGLE  DC  SQUID 


In  this  section  we  examine  the  effects  of  the  capacitance,  temperature,  and  resistance  of  two  junctions 
contained  in  a  ring  of  superconductive  material,  commonly  called  a  DC  Superconducting  Quantum  Inter¬ 
ference  Device  or  (SQUID).  Modeling  of  the  dynamics  of  the  single  DC  SQUID  and  DC  bi-SQUID,  as 
well  as  DC  SQUID  and  DC  bi-SQUID  arrays  in  both  one  and  two  dimensions,  without  temperature  noise 
or  the  capacitance  terms  can  be  found  in  [16]. 

3.1  SYSTEM  OF  EQUATIONS  FOR  A  SINGLE  SQUID 

Figure  10  is  the  schematic  diagram  of  a  single  DC  SQUID. 

u 
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Figure  10.  Single  DC  SQUID  schematic. 
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If  we  write  out  the  current 
tions: 

i  = 
h  = 

H  = 

Around  the  SQUID,  there  is  the  relation 

ifl  +  —ll  —  ^e~\~  ^2  +  ■^*2-  (2) 

To  get  an  equation  for  the  first  junction,  we  substitute  ^2  =  ~  U  into  Equation  (2): 

-ll  =  ^Pe  +  <^2  -  'fl  +  2  -  H)- 

Rearranging  gives 

R-  -L  -L 

PH  =  +  <Pe  +  P’2- 

h  —  ^  +  P>2  —  P>l)- 


and  phase  relationships  in  the  schematics,  we  get  the  following  six  equa- 


12  —  ^4,  is  =  n  +  ^2 

2eHiBT  ^ 


n  =  ^3, 
fo 

difi  . 

Pc +tci  sin  ipi 

ar^  dr  I  oh 

(fip2  d(p2  ,  .  .  2eKBT 


25{t  —  t) 
25  [t  —  t). 


1 


Substituting  ii  =  Z3  and  then  =  pc  .  o  H — ; - h  id  sin(^i - - — 2o[t  —  r  )  yields 

dr^  dr  I^Ti 


Pc 


d'^ifi  dip 

dr‘^  dr 

d‘^ipi  difi  iB  1 


a-ifi  aipi  .  2eKBT  ,  ,1  /  ,  x 

PC-jz^  +  —  +  ^ci  sin(^i - -  r  )  =  —  +  -((^e  +  ^2  -  ^i) 


2  /3 

2eKjBT 


dr^  +^-y  +  ^(^e  +  ¥>2-^i)-icisin^i+ 


25(t  —  r). 


Similarly,  we  can  substitute  i\  —  iB  —  ^2  into  Equation  (2)  to  find  a  second  governing  equation 


-  h)  =  +  ^2  -  ^1  +  ^^2. 


Rearranging  gives 


-Ph  =  +  V^e  +  ^2  -  ^1 

*2  =  y  -  ^(<^e  +  <^2  -y^l). 

Substituting  ^2  =  ia  and  then  +  *c2  sin (^2  —  25 (t  —  t')  yields 

dr^  dr  I^ri 

d?ip2  d(f2..  .  2eKBT  ,  is  1 

l3c^  +  ^  +  ^.2Sm.e2-^^2S(r-r)  =  ---(^,  +  ^,-.p,) 

d‘^if2  .  dip2  is  I  ,  X  .  .  2eKBT  , 

where  ipi  and  ip2  are  the  phases  across  each  of  the  Josephson  junctions,  (3  is  the  normalized  inductance  and 

Be 

ipe  =  2raxe.  The  parameter  Xe  =  —  is  the  normalized  external  magnetic  flux  per  unit  area,  and  a  is  the 

$0 

SQUID  area.  We  use  the  approximate  assumption  that  a  =  /3. 

These  two  second-order  differential  equations  can  be  written  as  a  system  of  four  first-order  ODEs: 
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26{r  —  r). 
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3.2  IMAGES  FOR  A  SINGLE  SQUID 

We  integrate  the  system  of  equation  in  Equation  (3),  varying  the  bias  current  over  time,  then  varying 
the  external  magnetic  field  over  time.  Eigure  1 1  shows  the  resulting  four  plots  with  an  ideal  response. 
Here,  Kp  =  0.01,  R  =  1.0,  C  =  0.1,  Xe  =  0,  and  (3  =  1.0  for  the  top  two  images,  and  Kp  =  0.01, 
R  —  1.0,  C  —  0.1,  ip  —  2.0,  and  [3  —  1.0  for  the  bottom  two  images.  The  top  two  images  are  basically 
the  same  as  the  first  case  for  the  single  junction,  but  now  there  are  two  junctions,  we  vary  ip  from  —10.0 
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to  10.0  and  back  to  —10.0.  The  bottom  two  images  are  the  voltage  over  time  and  the  average  voltage  re¬ 
sponse  when  sweeping  over  the  external  magnetic  field  from  —5.0  to  5.0  and  back  to  —5.0. 


Figure  11.  curve  (top  left)  and  I-V  curve  (top  right)  with  Kp  =  0.01,  R  =  1.0,  C  = 

Xe  =  0,  and  p  =  1.0.  V{xeR)  curve  (bottom  left)  and  V{xe)  curve  (bottom  right)  with  Kp  = 
R  =  1.0,  C  =  0.1,  ip  =  2.0,  and  p  =  1.0. 


0.1, 

0.01, 


For  the  second  simulation  presented  in  Figure  12,  we  increased  the  capacitance  to  (7  =  2.0,  while 
keeping  the  other  values  the  same  (Ffp  =  0.01,  i?  =  1.0,  Xg  =  0,  and /3  =  1.0  for  the  top  two  images, 
and  Kp  —  0.01,  R  —  1.0,  ib  —  2.0,  and  /3  =  1.0  for  the  bottom).  The  hysteretic  behavior  is  now  present 
in  the  /-l/curve.  The  presence  of  hysteresis  completely  destroys  the  SQUID  average  voltage  response. 
While  it  is  still  oscillatory,  it  looses  all  of  its  range  in  voltage.  Instead  of  stretching  from  0  to  0.9,  it  only 
oscillates  between  0.8  and  0.9.  This  negatively  impacts  sensing  signals  since  we  want  the  largest  voltage 
swing  possible. 


Figure  12.  V{IbR)  curve  (top  left)  and  I-V  curve  (top  right)  with  Kp  =  0.01,  R  =  1.0,  C 
Xe  =  0,  and  p  =  1.0.  V{xeR)  curve  (bottom  left)  and  V{xe)  curve  (bottom  right)  with  Kp 
R  =  1.0,  c  =  2.0,  ip  =  2.0,  and  p  =  1.0. 


2.0, 

0.01, 
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If  we  decrease  the  bias  current  to  be  closer  to  the  critical  current  of  the  inner  part  of  the  hysteric  loop 
on  the  I-V curve  (is  ^  1.55)  then  we  again  get  a  larger  range  from  0  to  0.65,  but  it  is  far  from  ideal  and 
there  is  a  large  amount  of  hysteresis  present.  The  other  parameters  in  Figure  13  are  iTp  =  0.01,  i?  =  1.0, 
and  C  =  2.0. 


Figure  13.  V{xe) 
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curve  with  Kp  =  0.01,  R  =  1.0,  C  =  2.0,  ip  =  1.55,  and  ^  =  1.0. 


For  an  even  larger  value  of  the  capacitance  (C  =  20),  the  average  voltage  response  is  non-oscillatory 
and  is  a  constant  line  at  V (xg)  =  1  (see  Figure  14). 


Figure  14.  V{IbR)  curve  (top  left)  and  I-V  curve  (top  right)  with  Kp  =  0.01,  R  =  1.0,  C 
Xe  =  0,  and  p  =  1.0.  V{xeR)  curve  (bottom  left)  and  V{xe)  curve  (bottom  right)  with  Kp 
R  =  1.0,  c  =  20,  ip  =  2.0,  and  P  =  1.0. 


20, 

0.01, 


If  we  go  back  to  capacitance  C  =  2.0,  while  decreasing  the  resistance  (i?  =  0.5)  with  the  other  values 
remaining  the  same  (Kp  =  0.01,  Xg  =  0,  and  /3  =  1.0  for  the  top  two  images,  and  Kp  —  0.01,  ip  —  2.0, 
and  /3  =  1.0  for  the  bottom),  then  we  restore  the  original  average  voltage  response  (see  Figure  15).  Some 
damping  is  still  seen  in  the  voltage  over  time  response  for  the  I-V  curve  simulation;  however,  it  doesn’t 
seem  to  impact  the  average  voltage  response  curve.  This  is  good  because  it  means  if  there  is  an  unavoid¬ 
able  amount  of  hysteresis  caused  by  a  capacitance,  then  the  response  can  be  improved  by  decreasing  the 
junction  normal  resistance.  Since  R  (x  C^,  while  we  simulated  the  response  we  didn’t  include  any  results 
here  because  it  looked  very  similar  to  Figure  14.  The  average  voltage  response  was  completely  destroyed 
and  there  was  hysteresis  in  the  I-V  curve. 
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Figure  15.  V{lB,t)  curve  (top  left)  and  I-V  curve  (top  right)  with  Kp  =  0.01,  R  =  0.5,  C  =  2.0, 

Xe  =  0,  and  /3  =  1.0.  V{xeR)  curve  (bottom  left)  and  V{xe)  curve  (bottom  right)  with  Kp  =  0.01, 

R  =  0.5,  C  =  2.0,  ip  =  2.0,  and  P  =  1.0. 

Next,  to  see  if  increasing  the  noise  temperature  improves  the  average  voltage  response,  we  set  Kp  = 
1.0  and  returned  to  the  case  where  C  =  2.0,  i?  =  1.0,  Xg  =  0,  and /3  =  1.0  for  the  top  two  images  of 
Figure  16,  and  C  —  2.0,  R  —  1.0,  ip  —  2.0,  and  /3  =  1.0  for  the  bottom.  The  hysteresis  is  reduced  slightly 
in  the  /-l/curve,  but  improvement  in  the  average  voltage  response  was  minor.  With  these  simulations,  we 
conclude  that  capacitance  is  bad  at  almost  any  amount,  and  care  should  be  taken  to  ensure  there  is  as  little 
as  possible  present. 


Figure  16.  curve  (top  left)  and  I-V  curve  (top  right)  with  Kp  =  1.0,  R  =  1.0,  C  =  2.0, 

Xe  =  0,  and  /?  =  1.0.  V{xe,t)  curve  (bottom  left)  and  V{xe)  curve  (bottom  right)  with  Kp  =  1.0, 
R  =  1.0,  C  =  2.0,  ip  =  2.0,  and  /?  =  1.0. 


To  see  the  effects  of  the  noise  temperature  on  a  desirable  average  voltage  response,  we  set  iTp  =  1.0 
with  C  —  0.1,  R  —  1.0,  Xe  —  0,  and  /?  =  1.0  for  the  top  two  images  in  Figure  17,  and  C  —  0.1,  Kp  —  1.0, 
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R  =  1.0,  iB  —  2.0,  and  (3  —  1.0  for  the  bottom.  There  is  not  much  change  in  the  I-V  behavior,  and  while 
the  V (xe)  curve  is  rougher,  the  overall  shape  was  maintained. 


Xe  =  0,  and  13  =  1.0.  V{xe,t)  curve  (bottom  left)  and  V{xe)  curve  (bottom  right)  with  Kp  =  1.0, 
R  =  1.0,  C  =  0.1,  iB  =  2.0,  and  p  =  1.0. 


The  last  simulation  we  performed  for  a  single  SQUID  is  a  case  that  could  be  considered  realistic  (see 
Figure  18).  There  was  some  capacitance  (C  =  0.5),  the  resistance  had  been  reduced  slightly  (R  =  0.9), 
and  there  was  some  noise  from  the  temperature  (Kp  =  0.2).  The  other  parameters  are  Xe  =  0,  P  =  1.0 
for  the  top  two  images,  and  ip  =  2.0,  P  —  1.0  for  the  bottom.  Here,  the  noise  is  low  enough  to  not  affect 
the  average  voltage  response  or  the  I-V  curve  much.  In  addition,  the  capacitance  is  offset  by  the  lowered 
resistance. 


Figure  18.  V(IbV)  curve  (top  left)  and  I-V  curve  (top  right)  with  Kp  =  1.0,  R  =  1.0,  C  =  0.1, 
Xe  =  0,  and  p  =  1.0.  V(xe,t)  curve  (bottom  left)  and  V(xe)  curve  (bottom  right)  with  Kp  =  1.0, 
R  =  1.0,  C  =  0.1,  iB  =  2.0,  and  p  =  1.0. 
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4.  SINGLE  DC  BI-SQUID 


For  our  study  we  are  interested  in  the  effects  of  having  two  ideal  junctions  around  the  loop  of  a  bi- 
SQUID  and  using  a  different  material  and/or  method  on  the  bisecting  Josephson  junction.  The  previous 
sections  have  explored  the  effects  of  temperature,  capacitance  and  resistance  on  the  basic  building  blocks 
of  a  bi-SQUID,  the  single  Josephson  junction  and  the  single  DC  SQUID.  In  this  section  we  will  explore 
the  single  DC  bi-SQUID  with  the  additional  terms  only  on  the  third  junction  in  the  same  manner.  First  we 
derive  the  phase  equations. 

4,1  DERIVATION  OF  EQUATIONS 

Figure  19  is  a  schematic  diagram  of  a  single  DC  bi-SQUID. 


I 

%  1 2 


Using  Kirchoff ’s  current  law,  we  get  the  following  six  current-phase  relations: 


— - Vicl  sin(/5i, 


d^2  ..  . 

^5  =  - |-^c2sm(^2 


il  +  ^2)  *1  +  *3 


=  *4)  h  =  *3  +  *5- 


Around  the  bi-SQUID  is  the  relation 


+  Llail  +  L2aH  —  +  -^16^2  +  -^2fe*5- 


(4) 


To  get  an  equation  for  the  first  junction,  we  substitute  i2  =  is  —  h  into  Equation  (4). 


+  Llah  +  L2aH  —  +  V^2  +  Llbi^B  “  U)  +  -^'26^5- 


Rearranging  gives 


{Lib  +  Lia)il  —  1/26*5  +  L2aH  —  +  ¥^2  “  +  Lihis- 
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Substituting  ii  =  —  iz  and  then  is  =  25{t  — 

dr^  dr  I^ri 

d(fi  ,  dip2  . 

— - h  ici  sin  (fi,  and  =  — - h  ic2  sm  (p2  yields 

dr  dr 


,  r  ^  •  /□  ^^3  d(f3  .  .  2enBT 

(Li5  +  Lia)  [  —  +  ^cl  sm  (fi  -  f3c  - -  ic2>  sm  (^3  + 


25{r  — 


r  ,  d‘-P2  .  .  \  T  f  ^  . 

^26  1  +  ^c2  sm  (^2  1  +  L2a  I  +  ^cl  Sm  (^1  )  =  (^e  +  ^2  “  +  ^16^5 


~  ^26 “^7  +  (^16  +  ^la  +  L2a)-^ - (^16  +  ^la) 


=  Lifcis  +  +  ¥^2  -  (/?1  -  (-^16  +  Lia  +  -^>2a)*cl  sin(^i 

(2eBBT  \ 

ic3sin(/33 - +-^26^c2sin(^2. 


Similarly,  we  can  substitute  ii  =  ib  —  ^2  into  Equation  (4)  to  find  a  second  governing  equation 


+  LiailB  —  ^2)  +  L2aH  =  +  ^2  +  ^16^2  +  ^26^5- 


Rearranging  gives 


L2aH  —  {Lia  +  Lii))i2  —  ^^26^5  =  +  ^2  “  ^1  —  LialB- 


Substituting  ^2  =  ^3  +  is  and  then  ^3  =  Pc—!^  + 

d,r^  d,r 


2eHiBT  ^ 


difi  .  .  dif2  .  . 

— - h^clsm(^l,^5  =  — - h^c2Sl 

dr  dr 


_  +  _  +  ,,3  3.„  ^3  -  (r  -  T 


sin  ip2  yields 

,T  ,  T  \(a  dip:i  .  2eKBT 

—  \Lia  +  Lii))  (  Pc  ^^2  ^c3  sm  ip3 


25(r  -  r')  +  ^  +  ic2sin(^2 

7o/i  ar 

ic2  sin  =P>e  +  <f2-<-Pl-  LlalB 


+  L2a  (  +  Icl  sm  Lpi\  -  L2b  I  +  t. 


L2a 


^  -  (Lia  +  Life  +  L2fe)^  -  (Lia  +  Life)  (/?c^  + 

=  —LialB  +  +  ^2  ~  ^1  +  (^la  +  ^16  +  ^26)^2  sin  (/^2 

/  2enBT  \ 

+  (Li«  +  Li5)  fic3sin(^3 - )  “  ^2a^ci  si: 


For  the  third  junction, 


^1  +  ^3  +  L2aH  +  (^3^  +  L^b)!^  —  iv26^5  +  ^2- 
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Substituting  ^3  = 


(f(p3 

dT^ 


— - 1-  ic2  sin  (f2  yields 

dr 


dips 

+  - h  ^c3  sm  (ps 

dr 


‘IcKjbT 

i^n 


26{r 


/ 


—  T 


,  i4 


dipt  ^  . 

— - h  ^cl  siny?i,  Z5  = 

dr 


(T  ,  r  «  ^^<^3  dips  .  2enBT 

{Lsa  +  Lsh)  [I5c^  +  —  +  t,s  sm  <^3  - 


25{t  — 


+  L2a  -J - h  Id  sm  ifi 

dr 


^26  1  +  ^c2  sm  (^2  )=  ^2  -  -  ^3 


{L^a  + 


dr  J  dr 

2eHiBT 


dr 


^2  -  -  ^3 


(2eKBT  \ 

ics  sin  ips - 2S(t  -t')]  -  L2aici  sin  ipi  +  L2bic2  sin  (p2 . 


Together,  we  have  a  system  of  three  differential  equations: 


-L2b^  +  (Life  +  Lia  +  L2a)^  -  (Life  +  Li„)  +  ^ 

=  Lii)iB  +  +  ^2  —  —  {Lib  +  Lia  +  i^2a)^cl  sin  (fi 

2eKjBL  ^ 


+ 


(2eKBT  \ 

ics  simps - )  + -^26^2  sin  (/32. 


L2a-^ - {Lia  +  Lib  + 

dr 


—  —LiaiB  +  +  ^2  ~  +  {Lia  +  Lib  +  L2b)ic2  Sm  (^2 

+  {Lia  +  Lib)  (^ic3  sin  cps  -  -  r  -  L2aid  sin  cpi 

fT  ^  T  L\  ( R  ^  dip^\  difl  d(f2 

[Lsa  +  Lsb)  {PC-^  +  -^  1  +  L2a-^  -^26-^  =^2-^1  -^3 

(2eK,BT  \ 

ic3  sin  (fs - 26{r  -r)j  -  L2aid  sin  (^i  +  L2bic2  sin  (p2 , 


where  the  normalized  critical  current  on  the  third  junction  is  ic3  = 
inductance  and  we  use  the  approximate  assumption  that  a  =  Lia  + 


,  L  =  27r^  is  the  normalized 

h  ^0 


Lib  +  L2a  +  ^26- 
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Splitting  the  second-order  differential  equation  gives  a  system  of  four  first-order  ODEs: 


—  Liijis  +  +  ^2  ~  ~  {LiJj  -\-  Lia  +  L2a)icl  sin  (pi 


—  —Liais  +  +  ^2  —  +  {Lla  +  ^16  +  ^26)^c2  sin  (/P2 


dips 


-  L2aici  sin  ifi  -h  L2bic2  sin  (/P2. 


(5) 


4.2  SIMULATION  RESULTS  FOR  A  SINGLE  BI-SQUID 

Figure  20  is  the  first  example  simulation  of  Equation  (5)  with  Kp  =  0.01,  i?  =  1.0,  C  =  0.1,  Xe  =  0, 
and  ic3  —  0.  The  additional  terms  of  iTp  and  f3c  are  only  on  the  third  junction  of  the  single  DC  bi-SQUID. 
We  set  the  inductances  to  Lia  =  Lib  —  0.27,  L2a  =  L2b  =  0.24,  and  L^a  —  L^b  —  0.3.  The  first  image 
is  the  voltage  over  time  as  the  bias  current  ip  is  varied  from  —10.0  to  10.0  and  back  to  —10.0  of  each  of 
the  four  differential  equations.  The  second  image  is  the  voltage  over  time  of  the  measured  output  of  the 
bi-SQUID  as  the  bias  current  ip  is  varied  over  the  same  range.  The  third  image  is  the  I-V  curve  for  the 
bi-SQUID,  which  is  determined  by  filtering  the  time-dependent  voltage  and  plotting  against  the  input  bias 
current. 


Figure  20.  Ui(/p,  t)  curve  for  each  Josephson  junction  i  =  1, 2,3  (left),  U(/p,  t)  curve  for  the 
bi-SQUID  (middle),  and  I-V  curve  (right)  with  Kp  =  0.01,  R  =  1.0,  C  =  0.1,  Xe  =  0,  and  ics  =  0. 

The  three  plots  in  Figure  21  are  the  average  voltage  responses  of  the  iCp  =  0.01,  i?  =  1.0,  C  =  0.1, 
ip  =  2.0,  and  ics  =  0.  The  first  image  are  the  voltage  over  time  plots  of  each  differential  equation  as 
the  external  field  Xe  is  varied  from  —10.0  to  10.0  and  back  to  —10.0.  The  middle  plot  is  the  voltage  vs. 
time  output  of  the  bi-SQUID  as  external  field  Xe  is  varied  from  —5.0  to  5.0  and  back  to  —5.0.  The  plot  on 
the  right  side  is  the  V  (xg)  curve,  which  is  determined  by  filtering  the  voltage  over  time  plot  and  plotting 
against  the  external  field. 


16 


Figure  21.  Vi{xe,t)  curve  for  each  Josephson  junction  i  =  1,2,3  (left),  V{xe,t)  curve  for  the 
bi-SQUID  (middle),  and  V{xe)  curve  (right)  with  Kp  =  0.01,  R  =  1.0,  C  =  0.1,  ip  =  2.0,  and  ic3  =  0. 


The  ideal  value  of  the  critical  current  on  the  third  junction,  ics  is  1.0.  The  average  voltage  response  for 
a  bi-SQUID  with  ics  =  1.0  is  shown  in  Figure  22.  This  is  the  value  for  which  V (xg)  is  the  most  linear. 


Xc 


Figure  22.  V(xe)  curve  (right)  with  Kp  =  0.01,  R  =  1.0,  C  =  0.1,  ip  =  2.0,  and  ics  =  1.0. 


Next,  we  increase  the  third-junction  critical  current  to  see  the  effects  in  the  I-V  and  V (xg)  curves  (see 
Figure  23).  The  top  three  plots  (here  Kp  =  0.01,  i?  =  1.0,  C  =  0.1,  Xe  =  3.0,  andics  =  0)  basically 
haven’t  changed.  But  the  V (xg)  plot  has  definitely  changed  since  hysteresis  is  present  (Kp  =  0.01,  R  = 
1.0,  C  —  0.1,  ip  —  2.0,  and  ics  =  3.0). 


Figure  23.  Vi(lB,t)  curve  for  each  junction  i  =  1,2,3  (top  left),  V{Ip,t)  curve  for  the  bi-SQUID  (top 
middle),  and  I-V  curve  (top  right)  with  Kp  =  0.01,  R  =  1.0,  C  =  0.1,  Xg  =  0,  and  ics  =  3.0.  Vi{xe,t) 
curves  (bottom  left),  V(xe,t)  curve  (bottom  middle),  and  V(xe)  curve  (bottom  right)  with  Kp  =  0.01, 
R  =  1.0,  C  =  0.1,  ip  =  2.0,  and  Zcs  =  3.0. 
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Increasing  the  capacitance  and  the  resistance  of  the  third  junction  to  (7  =  2.0  and  i?  =  2.0  yields  the 
six  plots  in  Figure  24.  All  other  parameters  stayed  the  same  as  the  prior  case.  While  the  amplitude  in  the 
top  left  side  image  and  the  spikes  in  bottom  left  side  image  of  the  red  and  cyan  plots  (relating  to  the  third 
junction  and  the  third  junction  second  derivative  differential  equations)  were  suppressed,  the  output  I-V 
and  V  (xe)  curves  were  not  affected  by  the  increase  in  capacitance  or  resistance  of  the  third  junction. 


Figure  24.  Vi{lB,t)  curve  for  each  junction  i  =  1, 2,3  (top  left),  curve  for  the  bi-SQUID  (top 

middle),  and  I-V  curve  (top  right)  with  Kp  =  0.01,  R  =  2.0,  C  =  2.0,  Xe  =  0,  and  ics  =  3.0.  Vi{xeR) 
curves  (bottom  left),  V{xeR)  curve  (bottom  middle),  and  V{xe)  curve  (bottom  right)  with  Kp  =  0.01, 
R  =  2.0,  C  =  2.0,  ip  =  2.0,  and  ^3  =  3.0. 


Increasing  the  capacitance  and  the  resistance  further  (to  (7  =  10  and  R  =  10)  results  in  a  smoothing 
of  the  I-V  curve  and  destroys  the  V (xg)  response  (see  Figure  25).  The  top  left  voltage  vs.  time  plot  has  a 
large  amplitude  on  the  cyan  curve,  which  is  transient  behavior. 


Time  xio*  Time  Xc 

Figure  25.  Vi{lB,t)  curve  for  each  junction  i  =  1,2,3  (top  left),  V{lB,t)  curve  for  the  bi-SQUID  (top 
middle),  and  I-V  curve  (top  right)  with  Kp  =  0.01,  R  =  10.0,  C  =  10.0,  Xe  =  0,  and  ic3  =  3.0. 
Vi{xe,t)  curves  (bottom  left),  V{xe,t)  curve  (bottom  middle),  and  V{xe)  curve  (bottom  right)  with 

Kp  =  0.01,  R  =  10.0,  C  =  10.0,  is  =  2.0,  and  ic3  =  3.0. 
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Returning  to  the  case  where  the  top  three  plots  in  Figure  26  have  =  1,67  =  0.1,  Xg  =  0,  and 
ic3  =  0,  and  the  bottom  three  plots  have  i?  =  1,  67  =  0.1,  =  2.0,  andigS  =  3.0,  but  increasing 

the  noise  temperature  parameter  {Kp  —  10)  has  almost  no  effect  on  the  I-V  curve.  The  amplitude  of  the 
red  curve  representing  the  voltage  over  time  of  the  third  junction  second  derivative  is  large  in  both  left  side 
images.  This  means  that  the  third  junction  can  be  a  significantly  noisier  junction  than  the  main  two  without 
affecting  the  I-V  curve.  The  V (xg)  curve  is  shaggier  and  looks  closer  to  an  average  voltage  response  with 
a  lower  ig3  value.  From  this  we  can  determine  if  the  method  we  are  using  to  create  the  third  junction  is 
noisy,  then  we  want  to  design  it  with  a  higher  ig3  than  we  would  if  it  was  a  theoretically  perfect  junction. 


Figure  26.  Vi{lBV)  curve  for  each  junction  i  =  1, 2,3  (top  left),  V{IbV)  curve  for  the  bi-SQUID  (top 
middle),  and  I-V  curve  (top  right)  with  Kp  =  lo,  =  1.0,  C  =  0.1,  Xe  =  0,  and  ^3  =  3.0.  Vi{xeV) 
curves  (bottom  left),  V{xeV)  curve  (bottom  middle),  and  V{xe)  curve  (bottom  right)  with  Kp  =  10, 
R  =  1.0,  67  =  0.1,  ip  =  2.0,  and  ics  =  3.0. 


Increasing  the  noise  temperature  parameter  further  {Kp  —  100)  makes  the  I-V  curve  simply  the  nor¬ 
mal  state  resistance  curve,  with  no  superconducting  area  (V  =  2/g),  shown  in  Figure  27.  The  V (xg)  curve 
no  longer  exhibits  oscillatory  behaviors.  This  is  most  likely  due  to  so  much  noise  on  the  third  junction 
(probably  an  unrealistic  amount)  that  it  is  leaking  into  the  other  junctions  and  driving  the  whole  bi-SQUID 
normal.  There  is  likely  a  value  of  Kp  between  10  and  100  for  which  the  third  junction  is  driven  normal 
(Zg3  to  zero),  and  the  dynamics  of  an  array  of  DC  SQUIDs  is  recovered.  More  investigation  is  needed. 

Finally,  we  look  at  the  effects  of  increasing  capacitance  and  resistance  in  the  presence  of  thermal  noise. 
The  top  three  images  in  Figure  28  have  Kp  —  5.0,  R  —  1.0,  67  =  0.1,  ip  —  2.0,  and  ig3  =  3.0,  and  the 
bottom  three  images  have  Kp  =  5.0,  R  =  2.0,  67  =  2.0,  ip  =  2.0,  and  Zg3  =  3.0.  The  increase  in  capac¬ 
itance  and  resistance  suppresses  the  noise.  Since  the  noise  is  suppressed,  the  effect  of  higher  linearity  with 
larger  ig3  values  is  reversed,  and  the  hysteresis  reappears. 
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Time  *10*  Time  *10*  Xc 

Figure  27.  Vi{lB,t)  curve  for  each  junction  i  =  1, 2,3  (top  left),  curve  for  the  bi-SQUID  (top 

middle),  and  I-V  curve  (top  right)  with  Kp  =  100,  R  =  1.0,  C  =  0.1,  Xe  =  0,  and  ic3  =  3.0.  Vi{xe,t) 
curves  (bottom  left),  V{xe,t)  curve  (bottom  middle),  and  V{xe)  curve  (bottom  right)  with  Kp  =  100, 
R  =  1.0,  C  =  0.1,  ip  =  2.0,  and  43  =  3.0. 


lime  xio  lime  xio*  Xc 

Figure  28.  Vi{xe,t)  curve  for  each  junction  i  =  1,2, 3  (top  left),  V{xe,t)  curve  for  the  bi-SQUID  (top 
middle),  and  V{xe)  curve  (top  right)  with  Kp  =  5.0,  R  =  1.0,  C  =  0.1,  ip  =  2.0,  and  ic3  =  3.0. 
Vi{xe,t)  curves  (bottom  left),  V{xe,t)  curve  (bottom  middle),  and  V{xe)  curve  (bottom  right)  with 

Kp  =  5.0,  R  =  2.0,  C  =  2.0,  ip  —  2.0,  and  *c3  —  3.0. 
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5.  SERIES  BI-SQUID  ARRAY 


Now  that  we  have  an  idea  of  how  the  single  DC  bi-SQUID  behaves  with  a  different  third  junction,  we 
explore  an  array  of  DC  bi-SQUIDs  coupled  together  in  a  series  fashion.  The  sensor  we  are  designing  will 
not  be  a  single  bi-SQUID,  but  an  array  of  thousands,  so  it  is  important  to  understand  the  effects  when  they 
are  coupled  together. 

5.1  SYSTEM  OF  EQUATIONS 


In  Figure  29  we  show  the  schematic  diagram  for  an  array  of  DC  bi-SQUIDs  coupled  in  series. 


Since  the  coupling  is  in  series,  we  don’t  need  to  derive  new  equations,  we  can  just  extend  the  single 
bi-SQUID  equations  to  an  array  of  N  bi-SQUIDs  by  including  a  mutual  coupling  term  of  strength  M  [16]. 
The  system  of  AN  first-order  ODEs  is 


{Llb^k  +  Lia^k  +  L2a,k)-^^  -  L2h,k-^^  “  {Llb^k  + 

~  Llb^k'^B  ^e,k  ^2^k  ^l,k  (-^16, “1“  ^la,k 

(2eK,BT  \ 

'^cS.k  sinc/93^^  vj^  —  25]^{t  —  t  )  J 

T  •  •  \ ^  f  C\  \ 

“1“  ^2bN^c2^k  Sm(/^2,fc  y  ^  ^^2,2  271X^0,1) 

i^k 

-  {Lia^k  +  Lib^k  +  L2b,k)-^^  -  {Lia^k  +  Lib^k) I^C.k-^ 

—  —Lia^k^B  +  ^e,k  +  ^2,k  ~  ^l,k  +  +  Lib^k  +  L2b,k)ic2,k  sin  (^2,fc 

(2eBBT  \ 

'^c3,k  sinc/93^^  vj^  —  25]^{t  —  t  )  J 

M 


^2a,k'^cl,k  SmC^X^^  -|-  y  ^  ^3^  V^2,z  ZTTXett^) 


i^k 


=  Vk 


d^3,k 

dr 

dipi^k  .  d>^2,k  ...  .  .  ..  dvk 

^2a,k—^ - ^2b,k—J^ - H  [Bsa^k  +  ^3b,k) PC,k-J^ 


~  ^2^k  ^l,k  ^3,k  dj2a^k'^c\^k  sinc^l^^  -|-  L2b^k'^c2^k  ^^^^2,k 

(2eBBT  \ 

ic3,k  sin  (^3^^  V]^  —  25^  (t  —  t  )  J  , 


(6) 


where  (pi^k^  ^2,k  and  are  the  phases  across  the  Josephson  junction  in  each  of  the  bi-SQUIDs,  k  = 

l,...,iV.  ’ 
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5.2  SERIES  BI-SQUID  SIMULATIONS 


To  start  off,  we  simulate  Equation  (6)  for  an  ideal  array  of  50  bi-SQUIDs.  This  array  has  small  capac¬ 
itance  (C  =  0.1),  small  temperature  noise  (Kp  =  0.01),  resistance  of  one  (i?  =  1.0),  and  third  junction 
current  of  one  (ic3  =  1)  (see  Figure  30).  We  set  the  midpoint  inductances  to  —  0.27, 

L2a,k  —  L2b^k  —  0.24,  L^a^k  =  ^36, =  0.3  and  vary  with  a  Gaussian  distribution  along  with  loop 
areas  between  a  =  0.52  and  a  =  1.52,  and  use  only  nearest  neighbor  coupling  with  M  =  0.001.  The 
top  left  image  is  the  voltage  over  time  of  the  array  as  the  bias  current  ip  is  varied  from  —10.0  to  10.0  and 
back  to  —10.0  when  Xg  =  0.  The  top  right  image  is  iho  I -V  curve  for  the  array,  which  is  determined  by 
filtering  the  time  dependent  voltage  and  plotting  against  the  input  bias  current.  The  bottom  left  image  is 
the  voltage  vs.  time  output  of  the  array  as  the  external  field  Xg  is  varied  from  —5.0  to  5.0  and  back  to  —5.0 
when  ip  =  2.0.  The  plot  on  the  bottom  right  is  the  V (xg)  curve  for  the  array,  which  is  determined  by  fil¬ 
tering  the  voltage  over  time  plot  and  plotting  against  the  external  field.  The  feature  we  are  interested  in  is 
the  anti-peak  around  the  zero  magnetic  field.  We  want  the  anti-peak  to  have  a  large  voltage  swing  and  be 
highly  linear  for  signal  detection  purposes. 


Figure  30.  E(/p,t)  curve  (top  left),  and  I-V  curve  (top  right)  with  Kp  =  0.01,  R  =  1.0,  C  =  0.1, 

=  1.0,  Xg  =  0,  =  50,  and  M  =  0.001.  V{xeR)  curve  (bottom  left),  and  V{xe)  curve  (bottom 

right)  with  Kp  =  0.01,  R  =  1.0,  C  =  0.1,  ic3  =  1-0,  ip  =  2.0,  N  =  50,  and  M  =  0.001.  The  loops  have 
a  Gaussian  distribution  between  a  =  0.52  and  a  =  1.52. 

Next  we  increase  the  capacitance  and  the  resistance  values  to  C  =  2.0  and  R  =  2.0  to  get  the  four 
plots  in  Figure  31.  All  other  parameters  stayed  the  same  as  the  prior  case.  The  output  I-V  and  V (xg) 
curves  were  not  affected  by  the  increase  in  capacitance  or  resistance  of  the  third  junction.  This  amount 
of  capacitance  and  resistance  on  the  third  junction  does  not  translate  to  hysteresis  in  the  I-V.  From  this  we 
can  determine  that  the  main  two  junctions  are  the  source  of  the  hysteric  behaviors. 

Increasing  the  capacitance  and  resistance  further  completely  destroys  the  anti-peak  as  shown  in  Figure 
32,  where  (7  =  10  and  R  =  10.  The  I-V  curve  still  shows  some  structure,  but  it  doesn’t  contain  a  segment 
where  there  is  zero  resistance  to  current.  There  is  still  not  much  hysteresis  present,  even  with  this  high  of  a 
capacitance  and  resistance  on  the  third  junction. 
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Figure  31.  curve  (top  left),  and  I-V  curve  (top  right)  with  Kp  =  0.01,  R  =  2.0,  C  =  2.0, 

ic3  =  1.0,  Xe  =  Q,  N  =  50,  and  M  =  0.001.  V{xe,t)  curve  (bottom  left),  and  V(xe)  curve  (bottom 
right)  with  Kp  =  0.01,  R  =  2.0,  C  =  2.0,  ic3  =  1.0,  ip  =  2.0,  N  =  50,  and  M  =  0.001.  The  loops  have 
a  Gaussian  distribution  between  a  =  0.52  and  a  =  1.52. 


Figure  32.  curve  (top  left),  and  I-V  curve  (top  right)  with  Kp  =  0.01,  R  =  10,C  =  10, 

ic3  =  1.0,  Xe  =  0,  N  =  50,  and  M  =  0.001.  V{xe,t)  curve  (bottom  left),  and  V(xe)  curve  (bottom 
right)  with  Kp  =  0.01,  i?  =  10,  C  =  10,  ^3  =  1-0,  ip  =  2.0,  N  =  50,  and  M  =  0.001.  The  loops  have 


a  Gaussian  distribution  between  a  =  0.52  and  a  =  1.52. 


Now  we  return  to  the  case  where  R  =  1.0,  and  C  =  0.1,  and  increase  Kp  to  10.  The  results  are 
in  Figure  33.  This  amount  of  noise  is  more  prevalent  in  the  voltage  over  time  plots  than  the  final  filtered 
response.  The  I-V  is  hardly  affected  at  all.  The  V [xg)  curve  has  some  loss  in  voltage  swing  from  5.0  to 
0.0. 
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Figure  33.  curve  (top  left),  and  I-V  curve  (top  right)  with  Kp  =  lo,  i?  =  1.0,  C  =  0.1, 

ic3  =  1.0,  Xe  =  0,  N  =  50,  and  M  =  0.001.  V{xe,t)  curve  (bottom  left),  and  V{xe)  curve  (bottom 
right)  with  Kp  =  10,7?  =  1.0,  C  =  0.1,  Zcs  =  1.0,  ip  =  2.0,  N  =  50,  and  M  =  0.001.  The  loops  have  a 
Gaussian  distribution  between  a  =  0.52  and  a  =  1.52. 


If  Kp  is  increased  to  100,  then  the  I-V  curve  is  driven  normal  and  the  anti-peak  disappears  com¬ 
pletely  (see  Figure  34).  This  is  the  response  we  would  see  if  the  bisecting  Josephson  junction  was  not 
superconducting,  while  the  main  two  junctions  were  superconducting.  The  bi-SQUID  array  does  not  just 
return  to  simple  SQUID  array  dynamics.  As  in  the  case  of  the  single  bi-SQUID,  this  is  probably  due  to 
the  noise  on  the  third  junction  measuring  so  incredibly  high  that  it  is  bleeding  into  the  other  junctions  and 
driving  the  whole  array  normal.  More  investigation  is  needed. 


Figure  34.  curve  (top  left),  and  I-V  curve  (top  right)  with  Kp  =  100,  R  =  1.0,  C  =  0.1, 

ic3  =  1.0,  Xe  =  0,  N  =  50,  and  M  =  0.001.  V(xe,t)  curve  (bottom  left),  and  V(xe)  curve  (bottom 
right)  with  Kp  =  100,  R  =  1.0,  C  =  0.1,  icz  =  1.0,  ip  =  2.0,  N  =  50,  and  M  =  0.001.  The  loops  have 
a  Gaussian  distribution  between  a  =  0.52  and  a  =  1.52. 
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The  final  simulation  we  ran  was  to  increase  all  three  parameters  a  little,  Kp  =  5.0,  i?  =  2.0,  and 
C  =  2.0  (see  Figure  35).  This  would  be  a  more  realistic  situation  with  a  poorly  constructed  third  junction 
and  near  perfect  main  junctions.  These  values  hardly  effect  the  anti-peak  of  the  V (xg)  curve  at  all.  The 
linearity  and  voltage  swing  seems  to  be  retained. 


Figure  35.  V{lB,t)  curve  (top  left)  and  I-V  curve  (top  right)  with  Kp  =  5.0,  R  =  2.0,  C  =  2.0, 
icz  =  1.0,  Xe  =  0,  N  =  50,  M  =  0.001.  V{xe,t)  curve  (bottom  left)  and  V{xe)  curve  (bottom  right) 
with  Kp  =  5.0,  R  =  2.0,  C  =  2.0,  icz  =  1.0,  ib  =  2.0,  N  =  50,  M  =  0.001.  The  loops  have  a 
Gaussian  distribution  between  a  =  0.52  and  a  =  1.52. 
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6.  FREQUENCY  DEPENDENCE  OF  A  BI-SQUID 


We  want  to  explore  the  effects  of  frequency  dependence  on  gain  and  linearity.  The  frequency  depen¬ 
dence  can  manifest  in  many  areas,  for  example,  the  conductivity  and  the  plasma  frequency.  The  conduc¬ 
tivity  is  the  inverse  of  resistance  and  the  plasma  frequency  depends  on  distances  between  the  junctions  and 
is  dependent  on  the  potential  carrier  density.  Here,  we  first  set  critical  current  of  the  third  junction  as  fre¬ 
quency  dependent  for  a  single  bi-SQUID,  with  schematic  as  seen  in  Figure  19. 


6.1  SIMPLIFIED  EQUATIONS 


Without  the  temperature  and  capacitance  terms,  the  single  bi-SQUIDs  equations  are 


dv?!  ^  d<^2  J  d(p3 

\-^l  “T  J^2a)  7  -^2b  7 

dr  dr  dr 


r  z'  r  r  ^  .  ^'^3 

^2a—; - +  ^2b)—; - 

dr  dr  dr 


¥^2  -  +  Luis  +  2nxea  +  Liics  sin  (^3 

+L26  sin  ip2  -  {Li  +  L2a)  sin  tpi 


ip2-<~Pi-  Liais  +  2TrXea  +  Liics  sin(/33 
-L2a  sin  ipi  +  {Li  +  L2h)  sin  Lp2 


difi 

^2a-r~ 

dr 


di^2 

dr 


r  d(p3 


^2  —  ^3  —  —  Lsics  sin  (^3 

-L2a  sin  (pi  +  L26  sin  (^2, 


where  Li  =  Lia  +  L16,  T3  =  L^a  +  ^36- 

The  external  signal  is  added  to  the  Xg  term  in  the  equations  as  A  sin(27r/r),  where  A  is  the  amplitude 
and  /  is  the  frequency.  To  simulate  the  frequency  dependence,  a  Heaviside  function  is  used  to  change  the 
value  of  ic3,  depending  of  the  incoming  signals  frequency.  Three  cases  were  explored.  One  was  that  the 
response  after  a  certain  frequency  would  act  like  a  traditional  SQUID,  or  ics  =  0  for  /  >  /g,  where  /g  is 
the  cut-off  frequency  for  the  Heaviside  function.  The  next  was  that  the  third  junction  line  would  act  like  a 
wire  connecting  the  two  sides  of  the  bi-SQUID.  We  represented  this  by  setting  ics  =  10.0  for  /  >  /g.  The 
last  case  explored  was  the  control  case,  where  there  was  no  frequency  dependence  introduced. 

6.2  FREQUENCY  SIMULATIONS 

The  average  voltage  response  of  a  single  SQUID  (blue)  and  single  bi-SQUID  (green)  is  shown  in  Fig¬ 
ure  36.  These  were  created  by  taking  the  average  of  the  time-dependent  voltage  response  at  a  point 

in  Xg.  This  procedure  was  done  for  a  range  of  Xg  with  =  2.002.  The  average  voltage  response  of  a  sin¬ 
gle  bi-SQUID  is  much  more  linear  than  an  single  SQUID.  This  is  important  for  signal  detection,  since  the 
greater  the  linearity,  the  smaller  the  spurious  signals. 

To  find  the  output  and  spurious  signal  strength,  the  time-dependent  voltage  was  determined  with  an 
input  signal  present,  then  a  power  spectral  density  was  taken  to  determine  the  output  signal  strength;  see  an 
example  at  300  MHz  in  Figure  37.  This  was  done  for  a  number  of  frequencies  in  the  30  MHz  and  3  GHz 
range.  For  each  of  the  three  cases  simulated,  the  normalized  critical  current  of  the  third  junction  was  equal 
to  I.O  until  the  cut-off  frequency. 
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Figure  36.  Average  voltage  response  of  a  single  SQUID  (blue)  and  a  single  bi-SQUID  (green). 


Figure  37.  Power  spectral  density  plot  for  ic3  =  1-0,  and  /  =  300  Mhz. 


The  output  signal  power  determined  from  the  simulations  of  a  frequency  input  between  30  MHz  and 
3  GHz  with  a  cut-off  frequency  Fc  =  500  MHz  for  a  signal  with  amplitude  0.25  is  shown  in  Figure  38. 
These  plots  were  created  with  Xe  =  0.25,  =  2.00.  The  control  case,  where  there  was  no  change  in 

ic3,  and  the  case  where  the  bi-SQUID  because  a  regular  SQUID  have  no  difference  is  output  signal  power; 
however,  the  case  where  the  critical  current  increases  to  10.0  the  output  signal  power  drops  by  more  than 
20  dB,  to  the  level  of  the  spurious  signals  in  Figure  39. 

The  largest  spurious  signals  can  be  seen  in  Figure  39.  The  largest  spur  increased  in  power  when  the 
critical  current  became  0.0  after  the  cut  off  frequency.  This  is  expected  due  to  the  difference  in  linearity 
in  the  V (xe)  between  the  SQUID  response  and  bi-SQUID  response.  For  the  case  when  the  critical  current 
increased  to  10.0  there  were  no  spurs  present  so  the  response  is  basically  at  the  noise  floor  of  the  system 
and  the  output  signal  power  itself  was  the  strength  of  a  spur  for  a  bi-SQUID  or  SQUID. 
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Figure  38.  Output  signal  power  vs.  frequency  for  a  input  signal  of  amplitude  0.25  using  a  Heaviside 
function  with  cut-off  frequency  fc  =  500  MHz  for  three  different  cases  when  /  >  fc\  control, 

ic3  =  0.0,  and  ic3  =  10.0. 


Figure  39.  Magnitudes  of  the  largest  spurs  vs.  frequency  for  a  Input  signal  of  amplitude  0.25  using  a 
Heaviside  function  with  cut-off  frequency  fc  =  500  MHz  for  three  different  cases  when  /  >  fc'. 
control,  ic3  =  0.0,  and  ic3  =  10.0. 


From  these  plots,  we  can  determine  that  a  frequency  dependence  can  be  modeled  and  visualized  for  a 
single  bi-SQUID.  The  next  step  is  to  increase  number  of  bi-SQUIDs  from  a  single  one  to  an  array. 
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7.  FREQUENCY  DEPENDENCE  OF  A  SERIES  BI-SQUID  ARRAY 


In  this  section,  we  expand  the  frequency  dependence  to  a  series  array  of  25  bi-SQUIDs,  and  this  time 
include  a  variance  in  the  time  scale  parameter  which  is  dependent  on  the  resistance  (or  inversely  the 
conductance). 


7.1  SERIES  ARRAY  MODEL 

The  system  of  equations  without  the  temperature  and  capacitance  terms  is 

/  7-  I  7-  \  Y  djip2^k  Y 

{Ll^k  +  L2a^k)^  - 

L/\h^k^B  ^e,k  ^l,k  {^l,k  ^2a,k)'^cl,k  sin(/9x^^  +  Ll^k'^c3,k 

-/IT  .  o  \ 

+  ^2b,k^c2,k^^^^2,k  +  2_^  ~  ~ 

i^k  ^  ^ 

~  Lla^k^B  ^e,k  “l“  ^2,k  ^l,k  {^l,k  ^2b,k)'^c2,k  ^^^y^2,k  ^l,k'^c3,k  sin(/93^^ 

AT  ,  c\  \ 

^2a,k'^cl,k  ^  ^  — 2  ^2^i  ZTlX^Cli) 

i^k 

dipik  j-  d(^2^k  Y  ^^3,k 

^2a,k  y  ^2b^k  V  r  J^3^k  V 

dr  dr  dr 

^2^k  ^l,k  y^3,k  ^2a,k'^cl,k  sin(/9x^fc  +  L2b,k'^c2,k  ^^^y^2,k  ^3,k'^c3,k  ^^^^3,k’) 

where  ^2,k  and  Lp^^k  are  the  phases  across  the  Josephson  junction  in  each  of  the  bi-SQUIDs,  k  — 

1, ...,  N  and  Li^k  =  Lia^k  +  L^^k  =  ^3a,)^c  +  L^b^k- 


7.2  SIMULATION  RESULTS 

First,  we  varied  the  frequency  and  held  ic3  and  ojq  constant  with  only  nearest  neighbor  coupling  for 
AT  =  0.001  (see  Figure  40).  Here,  Xe{f)  is  the  top  left  plot,  cJo{f)  is  the  top  right,  P{f)  is  the  bottom  left, 
and  AP(/)  is  the  bottom  right.  The  simulations  were  performed  with  Xe  =  0.25,  N  =  25,  and  AT  = 
0.001.  The  loops  have  a  Gaussian  distribution  between  a  =  0.52  and  a  =  1.52  and  the  frequency  is  in 
megahertz.  This  shows  some  differences  at  each  frequency,  even  if  all  other  parameters  are  the  same  and 
may  make  it  hard  to  determine  the  effects  of  each  element. 

To  see  if  there  was  any  distinct  structures  that  would  arise,  we  first  varied  ic3  and  ouq  together  (see  Fig¬ 
ure  41).  We  used  a  ramp-type  function  instead  of  the  Heaviside  function  from  the  last  section.  There  is  an 
overall  increase  in  power  of  the  signal,  and  an  overall  decrease  in  the  power  of  the  largest  spurious  signal. 
The  different  frequencies  still  vary  quite  a  bit,  however.  To  separate  the  variations  caused  by  frequency,  we 
held  the  frequency  constant  at  300  MHz  in  the  next  three  simulations. 

For  the  first  of  the  three  simulations  with  the  frequency  constant  at  300  Mhz,  we  varied  normalized 
parameter  uq  from  0.5  to  1000  and  held  ic3  constant  (see  Figure  42).  There  is  a  peak  around  uq  ^  1.0, 
then  a  decrease  until  ouq  ^  20  and  a  lower  peak  at  ouq  =  100.  The  largest  difference  in  power  between  the 
signal  and  the  largest  spur  is  at  ccq  ~  1-0. 
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Figure  40.  iczif)  (top  left),  uioif)  (top  right),  P(/)  (bottom  left),  and  AP(/)  (bottom  right)  with 
Xe  =  0.25,  N  =  25,  M  =  0.001,  and  a  Gaussian  distribution  between  a  =  0.52  and  a  =  1.52. 


Figure  41.  icsif)  (top  left),  wo(/)  (top  right),  P(/)  (bottom  left),  and  AP(/)  (bottom  right)  with 
Xe  =  0.25,  N  =  25,  M  =  0.001,  and  a  Gaussian  distribution  between  a  =  0.52  and  a  =  1.52. 


For  the  second  of  three  with  the  frequency  constant  at  300  Mhz,  we  varied  ics  from  0  to  1.0  and  held 
Wo  constant  (see  Figure  43).  The  signal  power  increases  and  largest  spur  power  decreases  continuously 
from  ic3  =  [0, 1.0],  resulting  in  the  largest  difference  at  ics  =  1.0.  Since  this  is  the  the  value  of  ics  for 
which  the  bi-SQUIDs  have  the  greatest  linearity,  this  was  expected. 


30 


Figure  42.  icsif)  (top  left),  Loo{f)  (top  right),  P(/)  (bottom  left),  and  AP(/)  (bottom  right)  with 
Xe  =  0.25,  N  =  25,  M  =  0.001,  and  a  Gaussian  distribution  between  a  =  0.52  and  a  =  1.52. 


Figure  43.  iczif)  (top  left),  uioif)  (top  right),  P(/)  (bottom  left),  and  AP(/)  (bottom  right)  with 
Xe  =  0.25,  N  =  25,  M  =  0.001,  and  a  Gaussian  distribution  between  a  =  0.52  and  a  =  1.52. 


Finally,  we  held  the  frequency  constant  at  300  Mhz  and  varied  both  ics  and  ujq  together  from  0  to  1.0 
and  from  0.5  to  1000,  respectively  (see  Figure  44).  The  ics  variations  dominated  the  outcome,  though  the 
notch  in  the  largest  spur  between  simulations  200  and  250  from  the  wq  variance  was  visible.  From  this  we 
can  conclude  that  the  value  of  the  third-junction  critical  current  has  the  largest  effect  on  the  signal  detec¬ 
tion,  and  this  is  what  should  be  focused  on,  control  through  the  frequency-dependent  methods. 
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Figure  44.  iczif)  (top  left),  wo(/)  (top  right),  P{f)  (bottom  left),  and  AP{f)  (bottom  right)  with 
Xe  =  0.25,  N  =  25,  M  =  0.001,  and  a  Gaussian  distribution  between  a  =  0.52  and  a  =  1.52. 


The  next  step  would  be  model  the  input  of  two  signals  at  once  and  to  increase  the  complexity  to  in¬ 
clude  the  capacitance  and  temperature  terms. 
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8.  EXPERIMENTAL  VALUES 


To  relate  the  model  to  realistic  values  that  can  be  experimentally  verified,  some  calculations  need  to  be 
performed. 


8.1  TIME  CONSTANT 


The  first  parameter  calculated  is  the  time  constant  that  normalizes  the  theory.  From  the  normalization 
of  the  phase  equations  that  model  a  single  bi-SQUID,  it  is  determined  that 

2e/oi?Ar 

T  =  (jjQl  =  - - - 1. 

h 


Values  for  each  of  the  parameters  can  be  gathered  or  calculated: 

Ic  ~  0.25mA  =  0.00025A  (example  value) 
Rn  ~  2AQ  (example  value) 

$0  =  =  2.067834  x  10“^^Wb 

2e 

h  =  6.6260696  x  10"^^Js 

=  1.0545717  x  10"^^Js 

27r 

e  =  1.6021766  x  10“^*^C 
Vc  =  IcRiv  =  eoo^tv  =  .0006V. 


Using  the  relations  and  values  above,  an  approximate  value  for  the  time  constant  can  be  determined: 


CCo 


2eI{)R]S[  2t{Vc 


27r.0006V 


=  1.8231208  X  10^^1/s  «  10^^1/s. 


V2i 


n 


$0  2.067834  x  lO-i^Wb 


Since  approximate  values  are  used  for  R  and  Ri^,  we  can  estimate  that  tuo  ~  10^^1/s.  This  means  that 
1  unit  of  r  is  equivalent  to  10“^^s  or  a  picosecond.  In  the  code,  values  of  r  are  used  from  0  to  2^®(ir  = 
1048.576,  since  dr  —  0.004.  Using  the  time  constant  determined  in  Equation  (7),  the  time  span  translates 
to  approximately  0  to  1.05  ns. 


8.2  OUTPUT  VOLTAGE  IN  VOLTS 

From  the  relations  in  Equation  (7),  the  response  of  the  bi-SQUID  in  volts  at  a  value  of  can  be  deter¬ 
mined  as 

V{t)  =  =  0.0003  {<pi  +  <p2)  Volts, 

since 

•  -^T/. 


8.3  AMPLITUDE  OF  INPUT  SIGNAL 

When  running  simulations,  the  minimum  detectable  signal  has  a  normalized  amplitude  of  A  =  0.000002, 
and  the  largest  signal  detected  without  distortion  has  A  =  0.25.  The  amplitude  is  in  Xe,  which  is  normal¬ 
ized  by  $0-  Thus,  we  can  calculate  the  flux  of  the  smallest  signal  detected  as  4.135668  x  10“^^ Wb  and  the 
flux  of  the  largest  signal  without  distortions  as  5.169585  x  10“^^Wb. 
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To  convert  this  value  to  volts,  we  determine  the  conversion  from  a  V (xe)  plot  of  a  single  bi-SQUID, 
which  has  no  gain.  For  Ic  =  0.00025A,  Rjy  =  2AQ  and  the  V (xg)  shown  in  Figure  36,  the  maximum 
voltage  is  14iax  =  0.0006  *  0.9V=  0.00054  volts.  The  V (xg)  goes  from  0  to  0.00054V  in  O.50o.  From 
this  we  can  determine  the  slope  of  the  V (xg)  as  5.223  x  lO^^VAVb.  To  get  the  signal  strengths  in  volts,  we 
multiply  their  flux  by  the  slope  of  the  single  bi-SQUID  AVR 

=  4.135668  X  10“^^Wb  *  5.223  x  lO^^VAVb  =  2.15  x  10“*^V  =  0.00215/iV 
=  5.169585  X  10“^®Wb  *  5.223  x  lO^WAVb  =  0.000269V  =  269/iV. 
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9.  CONCLUSION 


In  this  technical  report,  we  reviewed  the  effects  of  capacitance  and  noise  terms  for  a  single  junction 
and  a  single  SQUID.  Then  we  modeled  a  single  bi-SQUID  that  has  noise  temperature  and  capacitance 
terms  on  the  bisecting  junction  only.  Next,  we  expanded  to  an  array  of  50  bi-SQUIDs.  Finally,  we  look  at 
effects  of  a  frequency  dependence  of  the  third  junction  for  a  single  bi-SQUID  and  an  array  of  bi-SQUIDs. 
This  study  was  done  in  order  to  understand  the  dynamics  of  HTS  YBCO  bi-SQUIDs  and  to  investigate 
carefully  designing  the  junctions  to  control  the  linearity  and  voltage  swing  of  the  anti-peak.  The  linear¬ 
ity  and  voltage  swing  are  important  for  increased  gain  and  improved  accuracy,  and  frequency  dependence 
could  be  used  to  implement  filter-like  capabilities. 

From  Section  2,  we  determined  that  for  a  single  junction,  capacitance  is  detrimental  as  is  causes  hys¬ 
teresis.  Some  is  tolerable  if  there  is  also  a  decrease  in  the  junction  resistance,  since  this  counteracts  the 
hysteresis  arising  from  the  capacitance.  A  large  resistance  negatively  impacts  junction  performance  and 
has  an  even  stronger  impact  on  hysteresis  since  /3c  oc  while  /3c  oc  C.  A  small  amount  of  temperature 
noise  also  decreases  any  hysteresis  present;  however,  it  also  decreases  the  junction  critical  current  (Figure 
8).  Overall,  with  this  modeling,  for  the  single  junction  dynamics,  a  decreased  resistance  (or  alternatively 
an  increased  conductance)  is  desired. 

For  a  single  SQUID,  the  same  conclusions  as  for  a  single  junction  were  reached  in  Section  3.  Even  a 
small  amount  of  hysteresis  in  the  I-V  curve  destroyed  the  voltage  swing  in  the  V {x^)  curve,  which  is  the 
response  we  are  trying  to  manipulate  for  signal  detection.  Large  voltage  swing  equates  to  large  signal  gain 
and  linearity  equates  to  a  large  spurious  free  dynamic  range.  When  hysteresis  is  present  in  the  I-V  curve, 
the  bias  current  can  be  decreased  to  increase  the  voltage  swing;  however,  this  results  in  hysteresis  in  the 
V{xe)  curve  (Figure  13). 

When  we  just  include  the  capacitance  and  temperature  terms  on  the  bisecting  Josephson  junction,  and 
leave  the  outer  loop  main  junctions  as  ideal  junctions,  their  impact  on  the  V  (xg)  curve  are  decreased.  This 
was  explored  in  Section  4.  The  temperature  term  did  counteract  the  effects  of  a  large  third-junction  critical 
current  on  the  V (xg)  curve  (becomes  hysteretic)  when  there  was  small  capacitance  and  resistance  values 
(Figure  28).  When  the  capacitance  and  resistance  were  increased,  this  was  lessened.  The  effects  of  just 
modifying  the  third  junction  in  the  bi-SQUID  are  those  that  could  arise  from  junctions  either  fabricated 
using  a  different  method,  such  as  ion  beam  milling  or  ion  damage,  or  with  junctions  specifically  designed 
to  exploit  certain  properties,  as  with  multiferroic  junctions. 

In  Section  5,  we  explored  series  coupled  arrays  of  bi-SQUIDs  with  the  capacitance  and  temperature 
terms  on  the  bisecting  Josephson  junction,  and  leaving  the  outer  loop  main  junctions  as  ideal  junctions. 
Gain  is  increased  as  we  couple  the  SQUID  or  bi-SQUID  devices  together  into  arrays,  so  when  we  fabri¬ 
cate  the  prototype  bi-SQUID  sensor,  we  anticipate  having  tens  of  thousands  of  bi-SQUIDs  in  the  array.  As 
with  the  single  bi-SQUID  in  Section  4,  the  impact  of  the  capacitance  and  temperature  terms  on  the  V (xg) 
curve  was  decreased.  From  this  we  can  surmise  that  if  care  is  taken  with  the  main  two  Josephson  junctions, 
then  the  bisecting  junction  can  be  fabricated  with  an  inferior  method,  and  still  result  in  a  good  device.  An 
inferior  method  could  be  used  to  overcome  design  layout  issues  or  for  the  properties  gained  by  using  them 
(i.e.,  multiferroic s).  The  effects  on  the  linearity  of  the  anti-peak  still  need  to  be  explored. 

The  frequency  dependence  of  a  single  bi-SQUID  was  explored  in  Section  6.  Here  a  frequency  depen¬ 
dence  was  placed  on  the  third-junction  critical  current  in  the  form  of  a  Heaviside  function.  We  determined 
that  we  can  control  the  signal  detection  if  there  is  a  frequency  dependence  on  the  critical  current  of  the 
third  junction.  For  a  very  basic  example,  if  we  could  find  a  method  that  results  in  a  large  ig3  value  10) 
for  large  frequencies,  and  have  a  value  close  to  ig3  =  1.0  in  a  range  where  the  signal  strength  is  small,  and 
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^3  =  0  elsewhere,  we  could  have  a  device  that  acts  like  a  low-pass  filter  and  amplifies  the  signals  further 
in  the  determined  range. 

For  the  array  of  bi-SQUIDs  in  Section  7,  the  value  of  the  third-junction  critical  current  had  a  much 
larger  effect  on  the  frequency  dependence.  These  simulations  were  performed  with  different  arrays  in  each 
simulation.  To  get  a  better  comparison,  the  25  loop  sizes  should  be  held  constant  throughout  each  different 
run. 

This  study,  while  extensive,  is  by  no  means  exhaustive.  Many  configurations  and  parameter  variations 
either  were,  explored  but  not  presented  here,  or  have  not  been  attempted  yet. 
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